Splitting the voter criticality 
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Recently some two-dimensional models with double symmetric absorbing states were shown to 
share the same critical behaviour that was called the voter universality class. We show, that for an 
absorbing-states Potts model with finite but further than nearest neighbour range of interactions 
the critical point is splitted into two critical points: one of the Ising type, and the other of the 
directed percolation universality class. Similar splitting takes place in the three-dimensional nearest- 
neighbour model. 



Noncquilibrium phase transitions are recently attract- 
ing increasing theoretical interest. One of the motiva- 
tions driving the research in this field is the idea that, 
similarly to equilibrium systems, nonequilibrium contin- 
uous phase transitions can be also divided into several, 
hopefully not so many, universality classes [lj. Indeed, 
there are some examples that show that such a classifi- 
cation, at least to some extent, can be made. Particu- 
larly interesting results exist for models with absorbing 
states, and a notable example is the directed percola- 
tion (DP) universality class. A conjecture by Jenssen 
and Grassberger Q that models with a single absorb- 
ing state should belong to the same universality class 
has by now very convincing supports. Another univer- 
sality class encompasses models with double absorbing 
states that in one dimension (d = 1) includes also 

models which dynamics has some additional symmetry 
(parity conservation) |5j- However, it is also known that 
any asymmetry in the dynamics that would favor any of 
the absorbing states drives the system into DP universal- 
ity class [6j- Universality of models with more than two 
absorbing states is even more problematic since for mod- 
els with symmetric absorbing states some details of the 
dynamics might affect the critical behaviour too [j] . It is 
becoming clear that a task of classifyingnonequilibrium 
phase transitions is far from completed [8j . 

Numerical results that support the above classification 
come mainly from one dimensional models. The situation 
for d > 1 is less understood, but there are some results 
in this case too. In addition to numerous examples of 
models with a single absorbing state and thus belonging 
to DP universality class pj, there are some indications 
that for d > 1 models with more than two symmetric 
absorbing states the phase transition should be discon- 
tinuous 9]. In between, there are models with two ab- 
sorbing states. Dornic et al. shown that in this case and 
for d = 2 a group of models belongs to a new universal- 
ity class called a voter universality class ^} • Originally, 
a voter model was introduced in a rather nonphysical 
context of opinion spreading pd| . Later on numerous 
examples of related and physically more relevant mod- 
els were also studied 0] . Although the phase transition 
in the voter universality class is continuous, the decay 
of the order parameter p upon approaching the critical 



point is slower than any power law decay p ~ e", where 
e measures the distance from the critical point. Formally 
such a decay might be described with (3 = 0. In addi- 
tion, the time decay of the order parameter at criticality 
is also slower than any power law decay and is in fact 
logarithmic (p ~ l/ln(i)) as shown exactly Some 
other exponents of voter universality class were also de- 
termined fio| . 

An interesting feature of the voter model, that was not 
yet addressed, is the fact that at its critical point actually 
two phenomena seem to take place. One of them is the 
symmetry breaking between two competing states of the 
model, that is similar to the symmetry breaking in the 
equilibrium Ising model. The second phenomenon is the 
phase transition between active and absorbing phases of 
the model. Since there is no symmetry breaking tran- 
sition in one dimensional equilibrium Ising model, it is 
easy to understand that for models with two absorbing 
states in d — 1, the symmetry breaking must result from 
the absorbing phase transition and thus both phenom- 
ena take place simultaneously. But this is no longer the 
case for d > 1 models, and the coincidence of these two 
transitions should not be taken for granted. 

In the present paper we study a recently introduced 
nonequilibrium Potts model which dynamics has two ab- 
sorbing states . We show that the model with nearest- 
neighbour interactions on square lattice belongs to the 
voter universality class. However, with extended range 
of interactions (up to third nearest neighbours) the voter 
critical point is splitted. Starting from the disordered 
phase and reducing a temperature-like control parameter 
T, the model first undergoes symmetry breaking phase 
transition. Calculation of the Binder cumulant indicates 
that this transition belongs to the Ising universality class. 
Upon further decrease of T, the model undergoes a sec- 
ond phase transition into the absorbing phase. Since at 
this point the symmetry is already broken, this second 
transition, as expected, belongs to the DP universality 
class. Similar behaviour is observed for the nearest neigh- 
bour model in the three-dimensional case. 

That a certain type of a nonequilibrium critical point 
can be considered as superposition of two other critical 
points is in our opinion new, and hopefully, it should in- 
crease the understanding of nonequilibrium phase tran- 
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sitions. Let us notice that a superposition of different 
critical points exist in some equilibrium systems. For ex- 
ample, in the frustrated XY model two phase transitions 
of the Kosterlitz-Thouless type (U(l)) and of the Ising 
type (Z2) under certain conditions most likely happen 
simultaneously [T^j . Another example is a multicritical 
point in random magnets |l5j . 

Before defining our model let us recall that the equi- 
librium ferromagnetic 2-state Potts model can be defined 
using the following Hamiltonian ^(| 



H 



(hi) 



(1) 



where summation is over pairs of interacting sites 

on a cartesian lattice of the linear size L. With each site 
i we assign a variable = 0, 1 and S is the Kronecker 
delta function. 

To study model using Monte Carlo simulations 
one constructs a stochastic Markov process with suit- 
ably chosen transition rates. One of the possible choices 
corresponds to the so-called Metropolis algorithm. In 
this algorithm one looks at the energy difference AE be- 
tween the final and initial configuration and accept the 
move with probability min {l,e- AE / T }, where T is tem- 
perature. To obtain a model that would have symmetric 
absorbing states we modify the Metropolis algorithm of 
model Q as follows Q: when all neighbours of a given 
site i are in the same state as this site, then the site i can- 
not change its state. Thus, the dynamics of our model is 
defined as follows: (i) Select randomly the site i and its 
possible final state, (ii) If AE < z update the site i with 
the probability min{l, e~ AE ^ T }, where z is the number 
of neighbours interacting with site i (in our paper z is 
i-independent). Let us notice, that after such a modifi- 
cation, T is no longer temperature. Nevertheless, we will 
keep such a terminology. Moreover, the unit of time is 
defined as a single (on average) update of every site. 

To examine the properties of our model we used nu- 
merical simulations. A natural characteristic of models 
with absorbing states is the steady-state density of ac- 
tive sites p. In our model a given site i is active when 
at least one of its neighbours is in a state different than 
i. Otherwise the site i is called nonactive. Moreover, we 
used the so-called dynamic Monte Carlo method where 
one sets the system in the absorbing state with activity 
only locally initiated and measure some stochastic prop- 
erties of runs p"?! . Typical characteristics are the survival 
probability Pit) that activity survives at least until time 
t, and the average number of active sites N(t) (averaged 
over all runs). At criticality Pit) and N(t) are expected 
to have power-law decay: P(t) ~ t~ 5 and N(t) ~ t n . 
To detect a possible symmetry breaking in the model 
we measured the magnetization m = c^} — 1 and 

its variance s = Xd((2^i( cr i - — L d ) 2 ). At the 
symmetry breaking critical point the variance s, that is 
related with magnetic susceptibility for equilibrium sys- 
tems, is expected to diverge in the limit L — > 00. Below 
we present the results of our simulations. We ensured 




FIG. 1: The density of active sites p (+) and the inverse of 
the variance of magnetization s multiplied by a factor of 5 (x) 
as a function of temperature for the two dimensional nearest 
neighbour model. Simulations were done for L — 500. Close 
to the critical temperature T c (=1.7585 0) the linear decay of 



might be affected by some logarithmic corrections 
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that the lattice size L is large enough for the finite-size 
effects to be negligible. 

(A) Model with nearest-neighbour interactions on square 
lattice 

In this case the model was already studied by some of 
us, but only in the context of an absorbing phase tran- 
sition [9j|. A very slow decay of p at the critical point 
was observed, but being unaware of the relation with 
the voter model, we could not draw definite conclusions 
about the nature of the critical point in this case . In- 
deed, having a double symmetric absorbing state and be- 
ing driven only by interfacial noise, the model satisfies the 
Dornic et al.'s criterions for belonging to the voter univer- 
sality class. Additional confirmation is shown in Fig. ^ 
One can see that the variance s diverges at the same tem- 
perature where p vanishes. Moreover, s _1 seems to decay 
linearly at the critical point which implies that 7 = 1, as 
also predicted for this universality class |Toj . 

(B) Model with up to third-nearest-neighbour interac- 
tions on square lattice 

Studying our noncquilibrium Potts model within a mean- 
field approximation at the pair level, we noticed that for 
sufficiently large coordination number z a structure of 
the solution qualitatively changes in a way that clearly 
indicates two separate transitions in the model |19|. It 
prompted us to simulate our model with interactions in- 
cluding also further neighbours. First we studied the 
model with interactions up to second-nearest neighbours 
(z = 8). It turned out, that in this case either the model 
belongs to the voter universality class or there is only 
extremely small splitting beyond the resolution of our 
simulations. 

However, for the model with interactions up to third- 
nearest neighbours (z = 12), a qualitatively new picture 
emerges. Indeed, one can see in Fig.|5|that the variance s 
diverges at the temperature Tj that is clearly larger than 
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FIG. 2: The density of active sites p (+), the magnetization 
m (□), and the inverse of the variance of magnetization s 
multiplied by a factor of 30 (x) as a function of temperature 
for the two dimensional model with 12 neighbours (L — 400). 
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FIG. 3: The Binder cumulant U as a function of temperature 
T for the two dimensional model with 12 neighbours. The 
horizontal dotted line denotes the universal value U ~ 0.6107 
for the d = 2 Ising model. To diminish fluctuations long 
simulations were made with the simulation time t ~ 5 • 10 6 
Monte Carlo steps. 



the temperature T c where p vanishes. In the temperature 
interval T c < T < Ti our model is mag netized (Fig.©. 
To avoid rather slow coarsening effects in this interval it is 
better to start simulations from an asymmetric (m ^ 0) 
initial configuration. 

To examine the nature of the phase transition at 

T = Tf, we calculated the so-called Binder cumulant 

U = 1 — where m n is the n-th moment of magneti- 
sm^ ' & 

zation. In Fig. [3] one can see that at the crossing point 
U is relatively close to the universal value U — 0.6107 
of the d — 2 Ising model [2l]]. Although our model is 
nonequilibrium its critical behaviour at T — Tj is the 
same as in equilibrium (Ising) systems. Such a feature 
is in agreement with some expectations |22J that at the 
critical point of many nonequilibrium systems only some 
general properties (e.g., symmetry) determines the na- 
ture of the critical point while some other, as e.g., lack of 
detailed balance, are very often irrelevant in this respect. 

Since the symmetry is already broken upon approach- 
ing the critical point at T = T c , we expect that this 
critical point should belong to the (2+1) DP universality 
class. Simulations confirm these expectations. In Fig. 0] 
one can see that close to the critical point at T = T c the 
density of active sites behaves as p ~ (T — T C )P and we 
estimate that (3 = 0.61(4), that can be compared with 
the DP value 0.584 0. Additional confirmation is ob- 
tained using the dynamical Monte Carlo method, that 
enables us also to precisely locate the critical temper- 
ature T c — 4.7380(5). From the measurement of N(i) 
(Fig. |SJ at criticality we estimate rj — 0.25(3) which is 
in a reasonable agreement with DP value 0.230 [lj. We 
also measured P(t) and from these data (not presented 
here) we estimate 5 = 0.44(3) that can be compared with 
the DP value 0.451 [J. Our estimation of dynamical ex- 
ponents is much different from that obtained for some 
models of voter universality class S ~ 0.9, r\ ~ 0.0 

Let us notice that similarly to the nearest-neighbour 
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FIG. 4: The scaling of the density of active sites p in the 
vicinity of the critical point T c = 4.738 in the two dimensional 
model with 12 neighbours. Simulations were done for L = 500 
and the linear fit has a slope = 0.61. 
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FIG. 5: The average of number of active sites N(t) as a 
function of time t calculated using dynamical Monte Carlo for 
the two dimensional next-next-nearest neighbour model and 
(from top) T = 4.739, 4.7385, 4.738 (critical point), 4.7375, 
4.737, 4.736, 4.735. For each temperature the average is made 
usually over 2 ■ 10 4 independent runs (L = 5000). 
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FIG. 6: The density of active sites p (+), the magnetization 
m (□), and the inverse of the variance of magnetization s 
multiplied by a factor of 5 (x) as a function of temperature 
for the three dimensional nearest neighbour model (L = 60). 



case, the z = 12 model also has two absorbing states and 
is driven only by interfacial noise. Nevertheless, the criti- 
cal behaviour in this case is much different from the voter 
model. This is thus yet another example that shows that 
simple criterions most likely cannot be used for classify- 
ing noncquilibrium critical points. 

(C) Model with nearest-neighbour interactions on simple 
cubic lattice 

We also studied the three-dimensional nearest-neighbour 
version of our model (z=6). Similarly to the z = 12 case, 



here again the absorbing and magnetic phase transitions 
are separated (Fig. |SJ. We expect that in this case the 
magnetic phase transition belongs to the d = 3 Ising uni- 
versality class and the absorbing phase transition belongs 
to the (3+1) DP universality class. Confirmation of such 
a scenario will require, however, extensive numerical sim- 
ulations and is left for the future. 

In summary, we have shown that models with two ab- 
sorbing states in d > 1 dimensions might exhibit two 
transitions where the first one breaks the symmetry and 
the second one brings the model into an absorbing state. 
In the voter model and some related models these two 
transitions coincide. Hopefully, such an interpretation of 
the voter model will contribute to a better understanding 
of its unusual critical behaviour. For example, diminish- 
ing the strength of the further neighbour interactions we 
can reduce the splitting Tj — T c and examine a crossover 
to the voter universality class at which Tj = T c . Of 
course the Ising type phase transition is not the only type 
of the symmetry breaking and other types as e.g., 3-state 
Potts or XY, are also possible in nonequilibrium systems. 
It would be interesting to examine whether such critical 
points can be superposed with DP universality class, that 
might result in new types of critical behaviours. 
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